IMMUNE RESPONSE TO A MALARIA INFECTION: PROPERTIES OF A 

MATHEMATICAL MODEL 



Abstract. We establish some properties of a within host mathematical model of malaria 
proposed by Recker et al (6] [7] which includes the role of the immune system during the 
infection. The model accounts for the antigenic variation exhibited by the malaria parasite 
(P. falciparum). We show that the model can exhibit a wide variety of dynamical behav- 
iors. We provide criteria for global stability, competitive exclusion, and persistence. We also 
demonstrate that the disease equilibrium can be destabilized by non-symmetric cross-reactive 
responses. 
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pq \ 1. Introduction 

, This paper addresses the within host dynamics of the malaria parasite Plasmodium falciparum 

^ ' whose infection mechanism we briefly review here. Infection starts when a human is bitten by 

an infected mosquito that releases sporozoites in the bloodstream. The sporozoites quickly enter 
the liver where they mature, replicate, and differentiate into merozoites. The merozoites are 
O^, then released into the bloodstream, where they go on to infect erythrocytes (red blood cells). 

Merozoites reproduce within infected erythrocytes for a period of about two days. Finally the 
infected erythrocyte ruptures and releases new merozoites that repeat the infection cycle. A 
^ ' discussion leading to a mathematical model that considers a single parasite strain can be found 

, in [HI [3] and references therein. In practice however, there is a considerable diversity among the 

infected erythrocytes, which is reflected by a wide variety of the surface proteins (antigens) that 
are presented by the infected cells. A mathematical model that includes an arbitrary number 
of parasite strains was studied in a very elegant paper by Iggidr et al [S] where a competitive 
exclusion principle was established. Generically, only one strain survives while the others are 
driven to extinction. 



o 

, The mathematical models mentioned above do not include any immune response mounted by 

the human host. Although many details of immune responses to P. falciparum are presently 
not well understood, there is evidence that the antigenic variation between different strains of 
' the parasite prompts the immune system to mount both strain specific as well as cross reactive 

responses [U |6l [7] . The primary distinction between specific and cross- reactive responses is that 
they target major (unique to each strain) or minor (shared among strains) epitopes, respectively, 
on the infected cell's surface. 

The goal of this paper is to extend the analysis of the model proposed in [SI [7] which include 
the different immune responses described above. We provide some results concerning the global 
behavior of this model by 

(1) Showing global asymptotic stability of the system in two extreme cases (no cross immu- 
nity and perfect cross immunity). 

(2) Showing the possibility of oscillatory destabilization in the case of partial cross immunity. 

(3) Establishing conditions for both competitive exclusion as well as for persistence. 
Our results indicate that depending on parameter values, this model can exhibit a wide variety of 
dynamical behaviors. The full range of possible behaviors and biological implications is currently 
not fully understood and remains the objective of future research. 
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The rest of this paper is organized as foUows. In Section 2 we recall and slightly generalize 
the model from [6l [7] and in Section 3 we comment on the existence or non-existence of positive 
equilibria. In Section 4 we treat the case of a single parasitic strain and establish global asymp- 
totic stability, even when the growth rate of infected cells is assumed to be logistic, as opposed 
to linear. Similar results are obtained in Section 5 in two special cases: the case of no cross 
immunity, and the case of perfect cross immunity. In case of partial cross immunity, the dynamic 
picture is not as simple and this is illustrated in Section 6 by analyzing a particular example. 
In Section 7 we return to the general model and establish sufficient conditions for competitive 
exclusion as well as for persistence. These conditions are compared to similar ones for certain 
associated Lotka-Volterra systems of lower dimension. 



2. General modeling assumptions 

The model that we study here was originally proposed by Recker et al [6] and later analyzed 
by Recker and Gupta [7]. The model has the following form 

(2.1) yi = yi{(l)-azi-a'wi), 

(2.2) Zi = f3y^ - nzi, 

n 

(2.3) Wi = P'y^^Cjjyj - fi'wj, 

i=i 

where i ~ l,...,n. The variables yi, Zi, and Wi represent the abundance of the erythrocytes 
which are infected by the i-th parasite, and the magnitudes of the specific and cross-reactive 
immune response respectively. We assume that the immune responses are induced proportion- 
ally to the parasitic load at the rates /3 and j3' . The coefficients ^ and /i' model the life-span 
of the corresponding immune responses. The efficiency of both responses are given by a and 
a' . The coefficient represents the maximal growth rate of the parasite. We assume that all 
kinetic parameters are equal for all strains. Finally, we assume that each strain has a distinct 
major epitope, but two different strains may share common minor epitopes. In the model, we 
incorporate this assumption by introducing the non-negative cross-reactivity matrix C such that 
Cij > if the strains i and j share the same epitope and = otherwise. In the sequel we will 
refer to some special cases for which we introduce the following terminology: 

(1) We say that there is no cross immunity when C = I. 

(2) We say that there is perfect cross immunity when C = I'l, where 1 = (1 . . . 1) e R. 

(3) Otherwise we say that there is partial cross immunity. 

For mathematical convenience, we perform a simple rescaling of the original variables and rewrite 
system ((^ - ((^ as 

(2.4) yj = yi{l~ Zi-Wi), 

(2.5) Zi = yi- fiiZi, 

n 

(2.6) m = by^ Cijyj - H2Wi, 

and define 71 := — and 70 := — . 

In case of P. falciparum, there is a natural carrying capacity given by the number of avail- 
able erythrocytes which can be infected by the parasite. Setting aside the possible effects of 
erythropoesis, we can assume that such carrying capacity is constant and modify the model 
accordingly, 

y 

(2.7) ijt = yi{l ^ ^ ^0. 

(2.8) = yi- /iiZj, 

n 

(2.9) Wi = b^Cijyj - fi2'Wi. 
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3. The positive equilibrium 

Using the vector notation, we can express the equihbrium conditions of p.4p — (j2.6p as follows: 
z* = 7iy*, and w* = 71 Cy*. The positive equilibrium must then satisfy the condition 

(3.1) 7iy*+72Cy* - 1. 

In case of perfect cross- reactivity, where Cy = 1 for all there exists a positive solution of the 
form 

* - 1 
Vi^y^ \ , i=l,...,n, 

71 + ^^72 

which corresponds to a positive equilibrium. 

The positive equilibrium does not always exist. For instance, letting n = 3, 71 = 72 = 1 and 

/ 1 1+e 

(3.2) C= 1 + e 1 1 + e 

\ 1+e 1 

the solution of (|3.ip is given by (y^ , 2/2 j 2^3 ) = 1/2(2 - (1 + e)^) (1 - e, -2e, 1 - e), which is non- 
negative for e = 0, positive for small negative e, and neither for small positive e. 

4. Global stability in case n — 1 
In the simplest case n ~ \, the model 

(4.1) y = y(l-z-w), 

(4.2) z = y-fiiz, 

(4.3) w = by - fi2W, 
admits a unique positive equilibrium 

7i 72 

,71 +72' 71 +72' 71 +72 
which is globally stable. To see this, we rewrite (|4.1H4.3p as 

y = y{{z* ~ + ~w)), 

z = {y -y*) - f^i{z ~ z*), 

w = b{y - y*) - ^2{w - w*), 

and define 



iy\z\w*) 



V= —ds+ {s-z*)ds + - is^w*)ds. 

Jy* ^ J z* ^ J w* 

The function V clearly has a unique global minimum at {y^ ^ z* ^w^). In addition, 

F=(y-y*)((z*-z) + (ii;*-7«)) + (z-z*)((y-y*)-Mi(z-z*))+i(w;-u;*)(%-y*)-/i2(w-w*)) 
which simplifies to 

tV / *\2 A^2 / *n2 

V = —fii[z — z ) —{w — w ) . 

b 

Clearly, the equilibrium [y* , z* ,w*) is the only invariant set in {V = 0}. LaSalle's invariance 
principle then implies global stability of (?/*, z*, w*). 

Assuming a carrying capacity for the infected cells we have a different model 

(4.4) y = y(l - - z-w), 

(4.5) z = y-fiiz, 

(4.6) w = by — ^2W. 

It is easy to see that the modified model (|4.4H4.6p also admits a unique positive equilibrium 
{y* , z* ,w*). Using the same function V as before, we observe that 

V = -^{y y*f -^l^{z- z*f - ^{w - w*f. 
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We conclude again that the positive equihbrium is globally asymptotically stable. 

5. Global stability in case n > 1 

When there are two or more strains present, they can be antigenically distinct (no cross- 
reactivity), or antigenically similar (perfect cross-reactivity, see above), or there may be partial 
cross-reactivity. In this section, we prove global convergence for the first two cases. We also show 
that adding a carrying capacity does not alter the conclusions. 

5.1. Perfect cross- reactivity without carrying capacity. The equations are 

(5.1) iji = yi{l- Zi-Wi), 

(5.2) Zi = j/i - ^llZi, 

n 

(5.3) Wi = b^yj - fj,2Wi, 

for i = 1, n and they admit a unique positive equilibrium. We observe that for all i,j 

Wi - W.J = -/l2(Wi - Wj), 

hence all pairwise differences Wi — Wj decay exponentially to zero. To make this argument formal, 
using w — wi and Wj — w + Uj for j ^ 1 we rewrite equations (|5.1H5.3[) as 

(5.4) yi = yi(l-zi-w), yj ^ yj(l - Zj ~ {w + Uj)), j ^ 1, 

(5.5) = j/i-^iZj, 

n 

(5.6) w = b^yj - fi2W, 

(5.7) Uj = -H2U-j, j ^ 1. 

Clearly, the system (|5.4H5.7p is asymptotic to the limiting system 

(5.8) y^ = yi{l-Zi-w), 

(5.9) Zi = yi- /ii^j, 

n 

(5.10) w = b'^^yj ~ H2W- 
The Lyapunov function for (|5.8H5.1(J|) has the form 

b 



v=Ylyj + ]; / is-^*)ds- 



Indeed, after simplifications, we find that 



v = -^,,Y.{z,-z:f~^{w-w*f, 

i=l 

and then global asymptotic stability follows from Lasalle's invariance principle. 



5.2. Perfect cross-reactivity with carrying capacity. The equations are 

1 " 

(5.11) y^ = ?/,(l - — - - Wj), 

(5.12) z, = yi-fiiz^, 



(5.13) m = by^yj- H2Wi, 
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for i = 1, and they admit a unique positive equilibrium. Arguing as before, we consider the 
limiting system 



(5.14) 

(5.15) 
(5.16) 



w 



1 " 



w) 



n 



w, 



for which the Lyapunov function is 

n 

Indeed, after simplifications, 



Vi 



Vi 



ds+ {s-z*)ds) + 



(s — w*) ds. 



i=i 



Ml 



Z.i-Z^) -—{w-w ) , 



implying global asymptotic stability of the positive equilibrium. 

6. Analysis of a specific case with n = 2 and partial cross immunity 
In this section, we consider the dynamics of the system with n — 2 and C given by 

Notice that the dynamics of this system also arises when restricting the system with C given by 
(13. 2p and e = to the invariant set {j/i — 2/3, zi — z^, wi — w^}. The resulting equations have 
the following form 

Vi = 2/i(l - z., - Wi), i = 1, 2, 

Zi = Vi-fiiZi, i = l,2, 
wi = &(yi + ^2) - 
W2 = 6(2yi + 2/2) - M2W2. 
We re-introduce the coefficients 71 :— I///1 and 72 := 6//i2- The Jacobian of the system is given 

by 



/ 1- 


Zi — Wi 


-yi 


-yi 








\ 




1 


-Ml 
















b 





-M2 


b 



















I - Z2-W2 


-y2 


-2/2 













1 


-Ml 





V 


26 








h 





-M2 / 



This model admits at most four equilibria: 

(1) The zero equilibrium Eqq always exists and is always unstable since the Jacobian J{Eoo) 
(not shown) has eigenvalues Ai,2 = 1, A3, 4 = —Mi; ^5,6 = — M2- 

(2) The semitrivial equilibrium 



E 



10 



1 



71 



72 



71 +72 71 + 72 71 + 72 



,0,0, 



272 



71 + 72 



_ 71-72 



always exists. The Jacobian J{Eio) (not shown) has eigenvalues A4 — 
— /ii, Ag — —fJ-2, and Ai,2.3 are eigenvalues of the matrix 
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From the preceding stability analysis in Section 3, we already know that 5ft(Ai,2,3) < 0. 
Using the Routh-Hurwitz criterion, it is not difBcult to show that in fact 5R(Ai_2,3) < 0. 
Hence, the stability of Eio is determined by the sign of A4. Specifically, Eio is (locally) 
stable if 71 < 72, and unstable if 71 > 72. 
(3) The semitrivial equilibrium 

'0,0,. ' 



71 



71+72 ' 



A2 



' ' ' 71 + 72 ' 7i + 72 ' 71 + 72 ' 71 + 72 , 
always exists. The Jacobian J(-Eoi) (not shown) has eigenvalues Ai 
— /ii, A3 = — /U2, and A4,5,6 are eigenvalues of the submatrix 

-2/2 -2/2 

1 -Ml 
b -At2 

As we argued previously, 3?(A4,5,6) < 0. Since Ai > 0, Eqi is always unstable. 
(4) The nontrivial equilibrium En exists if and only if 71 > 72, i.e. precisely when Eio is 
unstable. The (2/1,2/2) coordinates of En are given by 

71 7i - 72 



(6.1) 



yi 



(71+72)^-27!' (71+ 72)2 _ 272 



2 • 



The common denominator is positive iff 71 > {V2 — 1)72, and the numerator of 2/2 is 
positive iff 71 > 72- The Jacobian at Eu is given by 



JiEn) 



/ 





-yi 


-2/1 








\ 




1 


-w 
















b 





-m 


b 






















-2/2 


-2/2 













1 


-Ml 





V 


26 








b 





-fJ'2 J 



As we showed previously, 

det J(£'ii) = 2/i2/2MiAi2((7i 



■72j 



272) > 0, 



thus J(-Eii) cannot have zero eigenvalues. It turns out, that in the special case /txi = 
/Lt2 = A*, all six eigenvalues of J {En) have strictly negative real parts: If /xi = /U2 = A*, 
the characteristic polynomial of J{En) has the following form: 

p(A) = (M + A)2 + C(l + b){yi + 2/2) + 2/12/2(1 + 26 - b^)^ , 

where ^ = A(m + A). Clearly, two roots are given by Ai,2 = —A*. The remaining four roots 
can be obtained by solving the quadratic equation in ^. We have 

1 + 26-62' 1 + 26- 62' 

hence 6 G [0, 1). Substituting the values of 2/1 and j/2, we have 



.2 I Mi + b){2-b) ^ ^i^{l-b) 



0. 



1 + 26-62 1 + 26-62 
The discriminant of this equation is 

o (l + &)^(2-6)^-4(l-6)(l + 26 - 62) 
^ (1 + 26 - 62)2 

Simplifying the numerator, we find that 

2 fc'(3-6)2 



(1 + 26 - 62)2 



Hence the roots are 



6 = -M, 6 = -M 



1-6 
1 + 26 - 62 ■ 
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The corresponding lambdas are solutions of 

1-6 
1 + 26- 62 



Ag 4 + ^A3^4 + = 0, Ag g + ^Ag^e + " ^, = 0. 



The positivity of coefficients in the above quadratics implies that 5R( A3. 4^5.6) < 0. 

6.1. Destabilizing the nontrivial equilibrium. In this section, we show that there exist a 
nonempty set of parameter combinations such that En is unstable. To do so, we fix the value 
6=1 and let fj,i = s, fj,2 = ce where c > and e is small. Recalculating the equilibrium values, 
we find 

c^e c(c — l)e 



2c - 1 ' " c2 + 2c - 1 
The positive equilibrium exists for all e > if and only if c > 1. The Jacobian of interest has the 
form (j6.ip with j/i, 2/2, /^i, given above. The characteristic polynomial of J{s) has the form 

piz,e) = s'^aoic) + e^ai{c){l + 0{e))z + e^a2ic){l + 0{e))z'^ 

+£2a3(c)(l + 0{e))z^ + ea4{c){l + 0{e))z'^ + £05(0)2^ + z^ 



where 



ao(c) 
ai(c) 
a2(c) 
03 (c) 
£14 (c) 



c^ic- 1) 
c2 + 2c - 1 ' 
4c^(e- 1) 
(c2 + 2c- 1)2' 

2c3(c- 1) 

(c2 +2C- 1)2' 

3c(2c-l)(c3 + 3c2 + c-l) 

(c2 +2C- 1)2 ' 

2c(2c- 1) 



c2 + 2c - 1 ' 
asic) = 2(c+l). 

Since p{z, 0) = z^, ./(O) has a zero eigenvalue of multiplicity 6. Now we expand the roots of p in 
powers of e. First, we evaluate p{ke°',e) and find that the leading terms are 

p{ke°',e) = e'^ao{c) + ke^+°'ai{c){l + 0{e)) + k^e^+^°'a2{c){l + 0{e)) 

+k'^e^+'^''a3{c){l + 0(e)) + fc^£i+^"a4(c)(l + 0(e)) + k^e^+^'^a^ic) + fcV", 

Now we construct the Newton diagram, that is, 

n{a) = min(4, 3 + a, 2 + 2q:, 2 + 3a, 1 + 4q!, 1 + 5a, 6a), 

which has two positive vertices at (1/2,3) and (1,4). Hence, the leading power of z is either 
a = 1/2 or a = 1. 

• Case a = 1 corresponds to z = fce + o(e). To determine the value of fc, we set the leading 
terms of p{ke,e) equal to zero and obtain the equation ao(c) + fcai(c) + k^a2{c) = 0. 
Simplifying this equation, we find that it is equivalent to 

^'^''"^^ (2fc2 + 4cfc + (c2 + 2c-l)) = 0. 



(c2 + 2c- 1)2 

Since c > 1, the roots are 

c- 1 

fci.2 = -C± ^ 

V2 

which are both strictly negative. 
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• Case a = 1/2 corresponds to z = re^/^ + + 0(e). Expanding 'pire^l'^ + fe,e), we find 
up to the two lowest orders of e that 

p(r£^/^ + le, E) = e^r^ (02(0) + ai{c)r^ + r*^ 

+re'''^{ai{c) + 2la2{c) + r'^as{c) + 'ir'^lai{c) + r'^ar,{c) + QrH^ . 

Setting the term equal to zero, we find that either r = (in which case we are back 
to the previous step) or that r satisfies the biquadratic equation 

a2(c) + a4(c)r^ = 0, 

which is equivalent to 

2^{c - 1) + 2c(2c - l)(c2 + 2c - ly + y + 2c - ifr^ = 0. 
The discriminant of this equation 

D = 4c2(c2 + 2c - 1)2 ((2c - 1)2 - 2c(c - 1)) = 4c'{c'' + 2c - 1)^ (c^ + (c - 1)^) 
is clearly positive, and both roots 

„2 



= ^2^|^(-(2c- 1) ± V(2c-l)2-2c(c-l)) 
are strictly negative. Hence, we have two pairs of pure imaginary values for r: 

ri,2 = ±i\ 



r3,4 = ±i} 





1) + V(2c-1)2 


- 2c(c - 


1)) 


c2 + 2c - 1 


'c((2c- 


1)-V(2c-1)2 


- 2c(c - 


1)) 



c2 + 2c - 1 

Substituting each pair into the e^/^ term and setting it equal to zero, we obtain the 
corresponding values of I: 

ai{c) + rl.^<i:i(c) + /■i.2"5(c) 



h = - 



12 = 



2a2(c) +4rf^2«4(c) +6rf_2 '' 

ai(c) + r| 4a3(c) + r| 405(0) 



2a2(c) + 4r|_4a4(c) + 6r|4 ' 

At this point, we have established the existence of six distinct branches of eigenvalues for small 
e > 0: 

zi = kis + o{e), 

Z2 = k2S + o{s), 

-23,4 = /i£ + ri,2e^^2 + o(e), 

-25,6 = hs + r3^4£^^^ + o{e). 

The first two eigenvalues are real and negative for small e > 0, so it remains to show that either 
h or I2 may bo positive for some values of c. 

The sign of the expression 2a2 + Aair^ + 6r'* can bo determined as follows. Consider a cubic 
polynomial f{x) = 2x{a2 + a^x + x?) which has three simple zeros at 2 < 4 < 0. Since 
J{x) > for a; > 0, we have that ./'(r^^a), /'(O) > 0, and /'(r| 4) < 0. Thus' 

f{rl^) = 2a2(c)+4r2_2a4(c)+6r4 2 >0, 

firl^) = 2a2(c)+4r2 4a4(c)+644<0. 

Since the denominators of l\ and I2 have opposite signs, it suffices to show that the numerators 
have the same sign. That would imply that one of is positive. We claim that the numerators 
of li and I2 are strictly positive for all sufficiently large c. Indeed, lets investigate the asymptotic 
behavior of the roots of the quadratics Qi{x) = ai(c) + 03(0)2; + a5{c)x'^ and Q2{x) = 02 (c) + 
a4{c)x + x"^. 
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• Equation Qi = is equivalent (after dividing through by 2c) to 
2c3(c-l) 3(c- l/2)(c3 + 3c2 + c- 1) 



(c2 + 2c^l)2 (c2 + 2c-l) 



2 



+ (1 + l/c)a;2 = 0. 



As c ^ oo, the roots of this equation converge to the roots of 2 + 3a: + a;^ = 0, that is, 
a; = — 2orx = — 1. This foUows from the continuity of roots. 
• Similarly, as c ^ oo, the roots oi Q2 — converge to the roots of 2 + 4a; + a;^ = 0, that 
is, X — —2 ± \/2. An equivalent statement is that 

hm rf 2 = -'2 - V2, hm rg 4 = -2 + \/2. 

C — >CC ' C — 'OO ' 

Since —2 — \/2 < —2 < — 1 < —2 + \/2 (i.e. the roots of Qi are located between the roots of 
Q2), we conclude that the numerators of li and I2 are strictly positive for all sufficiently large 
values of c. Since the denominator of h (respectively I2) is positive(respectively negative), . we 
conclude that li > and ^2 < for all sufficiently large c. (Numerically, this happens as long as 
c > 2.46.) We summarize the results of this section in the following Lemma. 

Lemma 1. Let 6=1, /ii = e, fi2 — cs, and 



C = 



1 1 

2 1 



then there exist e* > and c* > 1 such that for all < e < s* and c > c* , the Jacobian at the 
positive equilibrium En has two real negative eigenvalues, and two pairs of complex eigenvalues 
with positive and negative real parts respectively. In particular, the equilibrium En is locally 
unstable with two-dimensional unstable manifold. 

7. Results on boundedness of solutions, competitive exclusion and persistence 
7.1. Boundedness of solutions. Without loss of generality, consider the scaled model 

(7.1) Vt = yi{l-Zi-Wi), 

(7.2) Zi = yi- fiiZi, 

n 

(7.3) Wi = by^ Cjjyj - H2Wi, 

and suppose that b, fii, fi2 > and cu > for all i. 

Theorem 1 All nonnegative solutions of ( |7.i[ ) — ( |7. are ultimately uniformly bounded. 

Proof. Without loss of generality, we may consider only positive solutions , that is yi{t) , Zi{t) , Wi{t) > 
0. First, it is clear that since t/i < yi, we have yi{t) < yi(0)e*. Hence, all solutions are defined 
for t>0. Next, we introduce the quantities ai = yi/{zi + wi) > 0. It follows that 

, ^ yi{l - Zi~ Wi){zi + Wi) - yi{yi + 6I]"=i '^-ijHj ~ l^i^i ~ A^2Wi) 



Clearly, this implies that 

Using the fact that 

we obtain the inequality 



{zi + WiY 



Oil < Q!j(l - (1 + bcyi)ai H ). 

Zi + Wi 



1 < max(//i, ^2), Zi,Wi>\), 



on < ai{l + max(/^i, ^2) - (1 + bcii)ai). 



Hence, < as long as Ui > a* :— ^+™^^(^^^■^'■^) _ Consequently, ai{t) < ai := max(ai(0), a*) 
for all t > 0. Equivalently, we have that yi{t) < ai{zi(t) + Wi{t)), which implies that 

y^<y^{l-^), t>Q. 

ai 
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Therefore, yi{t) is bounded for alH > 0. Finally, we have that 



lim sup (i) < a*, limsupj/i(<) < a*, limsupzi(t) < — , limsupmi(i) < — o 

t — >oo t — *oo t — >QO f^l t — *oo ^2 

7.2. Competitive exclusion. Let -fi — 1/ fii and 72 — b/112, and define A = 71/ + 72C. 

Theorem 2. Suppose that the following condition holds: 
(7.4) 3r £ {1, ...,n} : Vx > 0, Ax>l^ {Ax)r > 1, 

then for any positive solution yi{t), Zi{t),Wi{t) > of ( |7. jp — ( | 7. i^P , we have linit^oo yrit) — 0. 

In (|7.4p . the vector inequalities correspond to the order induced by the standard cone _R". 
Proof. Let (/(t)) = j f{s)ds denote the time-average of the function f{t). Then for any 
positive solution, we have that 



Boundedness of solutions implies that 



= > 0, i-^oo, 

= ^ > 0, t-^oo, 

limsup(?/i/y^) = limsupl - - {wi{t)) < 0. 

t — *oo t — >oo 

Without loss of generality, there exists a convex compact set K C i?" such that y{t) E K for all 
t > 0. The convexity of K implies that {y{t)) E K for all t > 0. Let K' be the compact set 

K' = {xeK : Ax > 1}. 

By (|7.4p . compactness of and continuity, there exists e > such that {Ax.)r > 1 + e for 
all X E K' . Also by continuity, there exists 5 > such that (Ax)^ > 1 + e/2 for all x in the 
(5- neighbor hood of K' . 

Now we analyze the averages more carefully. Since 

n 

I {Z^) - 71 {y^) 1^0, I (W,;) - Cy (y^) | ^ 0, 

J = l 

we have that 

n 

limsupl - {zi{t)) - {wi{t)) = limsupl - -fi{yi{t)) - 72 ^ (t/j(t)) < 0, 

t — ^00 t — >00 . T 

that is, 

liminf(A(y(t))), > 1 
t — ^00 

for aU i = 1, n. It follows that there exists T > such that dist((y(t)), if') < (5 for aU t > T. 
Therefore, {A{y{t)))r > 1 + e/2 for all t > T. This in turn imphes that there exists T' > such 
that 

{yAt)/yr{t)) = 1 - (Z,(<)) - (Writ)) < -e/4, t > T', 

or equivalently, 

< yr(o)exp(-rf/4), t>'r'. 

This clearly implies that limt_+oo yrit) = 0. O 
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7.3. Partial persistence. Let 

(7.5) X ^ f{x,y) 

(7.6) y = 9ix,y) 

be a forward complete system on X x Y := M" x M™. Wc say that (|7.5p — (|7.6p is x-partially 
(strongly uniformly) persistent if there is some (5 > so that for all {x,y) G int(M") x int(M™) 
there holds that 

lim inf Xi{t) > S, i — 1, . . . ,n. 

t — >oo 

Inspired by the persistence result in 0] we have 

Theorem 3 Assume that dX xY is forward invariant for ( | 7. 5[ ) — (77S^ , and suppose K C XxY 
is a compact absorbing set (thus every forward solution of ( |7.5p — ^7.6^ eventually enters and 
remains in K). Let P : X x Y —t be continuously differentiate and the restriction of P to 
dX xY bcQ, and positive elsewhere. Assume that there is a continuous function ip : X xY —tM. 
so that 

P 

(7.7) —=il)onXxY \ {dX x Y) 

If for all {x,y) G dX x Y , there is some T > such that: 

(7.8) {i:{x{T),y{T)))>Q, 

then \7.5\ — \ 7.6^ is x-partially persistent. 

The proof can be found in [2] and is omitted here. 

Remark 7.1. Note that a result similar to Theorem 12.2.2 in but now for system \ 7.5\ — ( f 7. 6[ ) , 
remains valid. It states that Theorem 3 remains true if condition ( |7.^[ ) holds just for {x,y) which 
are uj limit points of orbits in dX x Y. The proof is exactly the same as m [4j. 

We will apply Theorem 3 to prove a persistence result for the malaria model (|7.ip — (|7.3p . 
which we re- write in a more compact form first: 

(7.9) X = diag(X)[l-(J„ I„)Y], 

(7.10) Y = -disig{p,)Y + BX, 

where ( ^ ) eWlx M^", 1 = (1 . . . 1)' e M", ^ = (//i . . . /^i ■ • -^2)' e M^" and 



B 



bC 



Note that 9M" x M^" is forward invariant, and that there is a compact absorbing set K in 
X M^ri by Theorem 1. Let 

A={In In)diag-\fi)B. 
We will show the following: Theorem 4 // there is some p E mf(M" ) so that 

(7.11) p'[l-AX]>0, 

for all X for which )BX^ equilibrium of ( |7.gp — ( |7.j0p in x M^", then 

system <\7.9\ i — ( |7. JOp is persistent. 

Proof. The proof proceeds in two steps. We will first show that system (j7.9p — (|7.10p is 
X-partially persistent using Theorem 3 and Remark 17.11 Then we will show that the system 
(frgj) - (fTTOl) is persistent. 

Step 1. Let us first establish X-partial persistence for (|7.9p — (|7.10p . Define the continuously 
diffcrcntiablc (perhaps by multiplying the vector p by a sufficiently large positive scalar) function 
P -.W^x M2j> ^ [0,00): 

P{X,Y)=^UXr, 
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which is on dW^ x M^r^ and positive elsewhere. Note that jLlI holds on x K^n^ ^ ^2,1-^ 
with 

^iX,Y) =p'[l-{In In)Y] 

We claim that for aU Z = {X, Y) e dW]^ x R'^'\ there is some T > such that: 

(V(z(r))) > 0, 

from which X-partial persistence will follow using Theorem 3. We will do this by induction on r, 
the number of non-zero components of X. If r = 0, then X{t) = for all < > 0, hence Y(t) — > 
as t ^ +00, so that uj{Z) = {0}. But since is an equilibrium point of ~ (|7.10p . (|7.1ip 
holds with X = 0, and therefore our claim follows from Remark 1 7. II Assume that the claim has 
been established for r = 1 , . . . , rn — 1 but that X has m non-zero components (of course, m < n). 
Denote the indices of these components by /, a proper subset of {0, 1, . . . , n}. There are two 
cases to consider: 

Case 1. The solution Z{t) converges to the boundary of the set D = {(X Y) ^ R" x 
R^" I X, ^ for aU iel}. Then oj{Z) is contained in part of the boundary of x R^" where 
at most m — 1 components of X are non-zero. The conclusion of our claim then follows from 
Remark 17. II and the induction hypothesis. 

Case 2. The solution Z{t) does not converge to the boundary of D. Then there is some e > 
and an increasing sequence ^ 00 so that Xi{t]^) > e for all k and all i € I. For i ^ I we 
have that Xi(t) = for all t > and thus in particular for all t — tk- Consider the (bounded) 
sequences of averages {X(tk)) and {Y{tk)), which we may assume -by passing to a subsequence 
if necessary- converge to limits X and Y with the property that Xi > ii i G I and Xi — 
otherwise. Integrating (I7.10p between and tk, dividing by tk and letting tk ^ 00 yields: 

(7.12) = -diagY + BX. 

Consider now the dynamics of the components Xi with i G / as described by ()7.9p . In particular, 
dividing by Xi, integrating between and tk, dividing by tk and letting — > 00, and using (|7.12p 
yields: 

0=1- {AX)„ i G /. 

Since Xi = for alH ^ / we see that ( -1^, ~ ) is an equilibrium of (17. 9P — (I7.10p . Finally 

Vdiag {p.)BXJ 

notice that as ifc — > 00, we have that: 

{^k{Z{tk))) -^p'[l~AX], 

which is positive by (|7.1ip . This establishes our claim. 

Step 2. In Step 1 we have shown that (|7.9p — (|7.10p is X-partially persistent, so that for all 
solutions starting in int(R" ) x int(R^") there is some 5 > Q such that 

liminfX(t) > 51, 

t—>oo 

where the above vector inequahty should be interpreted componentwise. Then (j7.1Qp imphes 
that for all large t, we have that 

Y > -diag(^)y + -Bl. 



This implies that: 



liminf y(t) > -di'Ag~\fi)Bl, 

t — ^CXD 2 



where the vector on the right-hand side has positive components, which establishes persistence 
of (17:^- (TTTUD . o 
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Figure 1. An example of system (|7.13p with n = 2. Nullclines are the dashed 
Unes. The hatched region represents The crosses represent the equihbria, 
the triangle Eq — Ei — E2 represents C, hence C n -D+ 7^ 0. 



7.4. Discussion. It is interesting to compare our competitive exclusion result (Theorem 2) and 
our persistence result (Theorem 4) obtained in the previous subsections to corresponding results 
for the following lower dimensional Lotka-Volterra system: 

(7.13) X = diag(X)[l-AX] 

For this system we can easily prove the following competitive exclusion result, using similar 
arguments as those in the proof of Theorem 2. 

Lemma 3 Suppose that ( |7.^P holds for system l\7.13\j . Then for any solution x{t) of l \7.13^ in 
intCMJ^), there holds that Xr{t) — > as i — > 00. 

For system ()7.13|) , there is the following persistence result [4] . 

Lemma 4 // there is some p G mf(R" ) such that l\7.11l holds for all X which are equilibria 
of ( |y. in dW^, then system l\7.1S^ is persistent. 

In other words, our conditions under which system (|7.9p ~ (|7.10p exhibits competitive exclusion 
(see Theorem 2), respectively persistence (see Theorem 4) holds, are the same as for the reduced 
order system (j7.13|) . 

Finally, we can interpret conditions (|7.4p and (|7.1ip geometrically, and will see that they are 
not mutually exclusive. This implies that there are examples of system (|7.9p — (|7.10p which don't 
fit our conditions for either competitive exclusion or persistence. 

In M", define the closed convex set 

D = {x eW'll- Ax < 0}. 

The boundary of D is given by those points x in D for which 1 — {Ax)i = for some i. In this 
case we say that constraint i is active for x. Condition (17. 4p says that there must be a constraint 
r which is never active in R" . 

Although a geometric interpretation of condition (|7.1ip is not immediately clear, it has been 
shown in [1] that (17. lip is equivalent to the following condition which does have a clear geometric 
meaning. 

(7.14) cnD+^d), 

where C is the convex hull of the set of equilibria of (|7.13p in 9R" and _D+ — DO R" . 

To see that the exclusion condition (|7.4p and (|7.1ip (or the equivalent (I7.14p ) are not mutually 
exclusive, consider a system (|7.13p with n = 2 with nullclines given in figure 17.41 Clearly neither 
condition (|7.4p nor condition (|7.14p hold. It is well-known that this is an example of a bistable 
Lotka-Volterra system. The equilibrium in int(R^) is a saddle and every solution in int(R^) not 
on the stable manifold of the interior equilibrium converges to either Ei or E2. 
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